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The effect of a magnetic field on Josephson current has been studied for a 
superconductor/normal- metal/superconductor (SNS) system, where N is a two-dimensional 
electron gas in a confining potential. It is found that the dependence of Josephson currents 
on the magnetic field are sensitive to the width of the normal metal. If the normal metal is 
wide and contains many channels (subbands), the current on a weak magnetic field shows a 
dependence similar to a Fraunhofer-pattern and, as the field gets strong, it shows another type 
of oscillatory dependence on the field resulting from the Aharonov-Bohm interference between 
the edge states. As the number of channels decreases (i.e. normal metal gets narrower), how- 
ever, the dependence in the region of the weak field deviates from a clear Fraunhofer pattern 
and the amplitude of the oscillatory dependence in the region of the strong field is reduced. 

KEYWORDS: SNS system, two-dimensional Normal metal, Fraunhofer pattern, Edge state, quantum wire 



§1. Introduction 

It is well known that in a superconductor /insulator /superconductor (SIS) system, Josephson 
current as a function of a magnetic field shows a Fraunhofer pattern.^ This remarkable phenomenon 
is caused by the interference effect of tunneling Cooper pairs. 

In the case of a SNS system, however, where N is a two-dimensional normal metal with a meso- 
scopic scale, the dependences of the Josephson current on the magnetic field are expected to be not 
so simple. For example, if the field is strong, the normal metal becomes a quantum Hall system and 
the edge states are formed. Ma and Zyuzin studied a SNS system under such a strong magnetic 
fielc&i where N is two-dimensional normal metal confined by a linear potential and the edge states 
of the lowest Landau level are filled, and two superconducting leads are connected to the normal 
metal by point contact junctions. They proposed that Josephson current flows through the edge 
states in the normal metal and as a function of the field, the current shows Aharanov-Bohm type 
oscillations, similar to oscillations of the physical quantities in quantum dots. 

* E-mail: ishikawa@watson.phys.s.u-tokyo. ac.jp 
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In this paper, we study a SNS system where the normal metal is sandwiched between bulk 
superconductors. Recently, such a structure has been actually constructed by Takayanagi et al% 
An external magnetic field penetrating the system can be arbitrary. The normal metal is assumed 
to be under the parabolic confining potential in the direction perpendicular to the junction resulting 
in the discrete channels. It turns out that the dependences of Josephson current are sensitive to 
the number of channels, N, being defined as the number of channels in the absence of a magnetic 
field. 

The organization of this paper is as follows. In §2, the model and a derivation of Josephson 
current are presented. In §3, we summarize the general features of the dependences of the current 
on a magnetic field and N. §4 and §5 are devoted to a system with large N under a weak field 
and a strong field, respectively. In §6, we investigate how these results in a system with large N is 
affected as N decreases. In §7 we summarize these theoretical results and discuss its realizability. 

§2. Model and Derivation of the Josephson current 

We consider a system as shown in Fig. 1, where two bulk superconductors are attached to a 
two-dimensional normal metal under the confining potential. We assume that a magnetic field is 
perpendicular to the two-dimensional normal metal and penetrates only the normal metal region 
uniformly. 

The Hamiltonian of the system is 

H = HSi + 7~tS2 + 7~(-N + T^Tunnel , (2.1) 

where Hsi(i = 1)2) and TCn are the BCS Hamiltonian of two superconductors and the two- 
dimensional normal metal, respectively. 

In the normal metal, the confining potential is assumed to be parabolic as&& 

V{y) = -mca^y 2 , (2.2) 

where m is the electron mass, and the width of the normal metal in the y-direction, W, is classically 
defined as 

1 2 fW\ 



where \x is the chemical potential. In the x-direction, electrons can move free along the length L. 

The vector potential A is taken as A = —Hy e x in the normal metal, and A = in the 
superconductors. This choice of gauge enables us to regard the macroscopic phase of Cooper 
pairs in two superconductors as spatially uniform. In the calculation in this paper, for the sake of 
simplicity, we neglect Zeeman energy whose effects are briefly discussed in §7. 

For the model described above, an eigenfunction in the normal metal is given by 
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where n = 0,1,2,... and k = 0, ±7r/L, ±2ir/L, ... and H n is the Hermite polynomial and 
2/o (k) = kl^/l 2 , Iq = y/h/muj , l c = y/h/muj c , uj = \Juq + uj 2 , w c = eH/mc being the cyclotron 
frequency. We denote this eigenstate as {n,k}. 
The eigenenergy E n k of the state {n, k} is given by 

E nk = huj [n + - H , (2.5) 

V 2 J 2M ' v ; 

where M = m(u! / loo) 2 . 

This spectrum is shown in Fig. 2. The Fermi wave number of each channel (subband, specified 
by n) is given by 



1 1 — hx \ n + — 




(2.6) 

where fi 2 K 2 F /2M = fj, ( i.e. K F = Wj2l\ ). The chemical potential \x is determined to conserve 
the number of electrons in the normal metal under a varying magnetic field. Hence, as the field 
increases, the number of occupied channels decreases. 

The states around {n, ±/cp(n)} contribute to transport phenomena at low temperatures. The 
characteristics of the eigenfunction, ^ n _|_fc j; ,( n )(r), are as follows.^ The center of the wave function 
in the y-direction, yo(±kF(n)), is given by 



uj 2 y /x 

The spatial extent of ip n ±k F (n){ r ) m the y-direction around yo(kF(n)) is roug hly l ox /n + 1/2. In 
the region of a weak magnetic field (o; c /a;o < 1), yo(kF(n)) is small compared with the classically 
defined edge W/2. Therefore the currents carried by the states {n,/ci?(n)} and {n, — ^(n)} are 
almost canceled all over the normal metal. As the field is increased {u c /loq > 1), the factor u> c /u> 
in eq. (2.7) approaches to 1. So yo(±kF(n)) is written as 

and spatial extent becomes small and converges to the Larmor radius Z c . Then the difference 
between yo(kF(n)) and yo{—kF{n)) becomes large, and yo(±kF(n)) for a small n almost coincides 
with the classically defined edge ±W/2. So the currents carried by {n, kpin)} and {n, —kp^n)} 
are spatially separated and this leads to the formation of the edge states, that is to say, the states 
carrying current along +x shift to one side of the normal metal, while states carrying current in the 
other direction shift to the other side of the normal metal. These influences of a magnetic field on 
the eigenfunction are in accordance with a classical viewpoint, since the Lorentz force (— e)v x B 
is opposite for electrons moving in opposite directions. 
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So far are discussed electronic states in the normal metal. The tunneling process between the 
superconductors and normal metal is described by H-Tunnei i n ec i- (2-1) as the following form, 

^Tunnel = 9^2 J dr ^ S ( X )^S1 a ( r ) $N a 0) + 6(x - L)^ a (r)$ S 2 a (r) + /l.cj , (2.9) 

where ipsia( r ) an d ^Na{ r ) are the field operators of two superconductors and the normal metal, 
respectively. 

Josephson current is given by the general formula as follows^ , 

2edfl , 

J = Tm- < 2 ' 10) 

where VL and 6 are the thermodynamic potential and the phase difference between the two super- 
conductors attached to the normal metal, respectively. The dependence of Q on 9 can be obtained 
by a perturbative expansion with respect to the pair amplitude of superconductors We evalu- 
ate Josephson current to the lowest order, which means considering the most simple process of a 
Cooper pair propagation as shown in Fig. 3. This approximation for the Josephson current is valid 
under the condition L > £ c , where £ c is a coherence length in the normal metal, representing a 
characteristic penetration length of a Cooper pair into the normal metal, and given by £ c = vp/2irT 
in a clean limit. We evaluate the pair amplitude in the local approximation as 

Aj(r -r',r- t') =g 2 < i>Si\{r,r)^ S i i(r' ' ,t') > 

- 9 2 <i>Sif(r,T)ipsii(r,T') > S(r - r') , (2.11) 

which means that two electrons of a Cooper pair in the superconductors tunnel into the normal 
metal at the same point. 
Based on these approximations, we can express Josephson current by the following form, 

J =j^(-i)T^2\A(uj m ,H)\ 2 J dyi J dy r j g N ^(r h r r ;uj m )g Nl {r u r r ;-io m )e ie 

- GN<t(r r ,rr,u} m )g N i(r r ,ri; -uj m )e~ ie j , (2.12) 

where T[ = (0, yf) and r r = [L, y T ) are the coordinates at the interfaces between the superconductor 
and the normal metal, and Qn<j is a thermal Green's function in the normal metal. In eq. (2.12), 
|A(a> m ,iT)| is defined as 

M<o m ,H)\=g. (2.13) 

where A(iJ) is the order parameter in the superconductor under a magnetic field H. With the 
eigenfunction in eq. (2.4), Qno can be written as 

g Na (r,r';co m ) = £ ^M^l ( U = E nk - p ) . (2.14) 

nk 1 L ° m ^ nk 
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By use of this Green's function, eq. (2.12) is transformed into 
J =f T E |AK,,#)| 2 Re[| d m J 



U! m >0 

= J c sm9, 



sm( 



(2.15) 



where we took the symmetries of the Green's function into account (see Appendix A). This result 
coincides with a well-known form of Josephson current proportional to the sine of the phase differ- 
ence, while the higher order correction in this perturbation, which is ignored here, gives rise to the 
higher harmonics of sin 9. 

J c defined in eq. (2.15) can be calculated as follows ( see Appendix B ), 



J c = i T Y. |AK„,^)| 2 ReE 
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2vr {(q + 01 - a 2 ) 2 + (61 + 6 2 ) 2 } {(q + ai + a 2 ) 2 + (61 - 6 2 ) 2 } 

where q = k\ — k2 (the relative wave number in the x-direction), k\ = a\+ibi, K2 = —02+162 (ai, bi > 
0) and a = (q + 2a x )L - 1/2 • tan -1 ^™/^ - fiu(ni + 1/2))]. 
Ki and K2 are defined as 



(2.16) 



2M 

h 2 k 2 



J- = fi — Huj (n\ + 5) + iw. 



lUJ m 



is defined by 



( Im[Ki],Im[K 2 ] > ). 



M = max{ni, n 2 } 
N = min{ni, n 2 } 



(2.17) 



(2.18) 



where is the Laguerre polynomial. 
For channels ni satisfying fi — hu>(rii + 1/2) > , a» and bi can be approximated as 



(2.19) 



The damping term exp(— 2biL) in eq. (2.16) can be written as exp(— (2m + 1) • L/£ c (ni)) where 
£c(ni) = hvF(ni)/2TrT. Here we consider the case where £,c(n\) < L for all occupied channels m, 
therefore in the summation over the thermal frequency io m (> 0), the term corresponding tom = 
becomes dominant. 
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§3. General Features of the Dependences on Magnetic Field and N 

Fig. 4 indicates the general features of the dependence on the magnetic field for various N, the 
number of channels at H = 0, which is related with the classically estimated width W introduced 
in eq. (2.3) as follows 



Therefore the large N corresponds to the wide normal metal. As the field is increased, the number 
of channels decreases and finally becomes 1 as seen from Fig. 4. 

The dependence is roughly classified into two regions; " Fraunhofer pattern" region and " Edge 
current " region. In " Fraunhofer pattern " region, Josephson current shows a dependence on the 
field like a Fraunhofer pattern as in a SIS system and in " Edge current " region, the dependence 
shows another type of oscillation due to the interference between the edge states in the normal 
metal. 

§4. Dependence in the Region of Weak Magnetic Field and Large N 

In this section, we assume that N is large enough (N ~ 100) and consider the case of the 
weak field {uj c /ujq < 1). Since a weak magnetic field is expected not to affect so much an orbital 
motion of electrons in the normal metal, Josephson current will show a dependence similar to a 
Fraunhofer-pattern in a SIS system. 

The numerical result is shown in Fig. 5. In this calculation, parameters are chosen as W = 
28 £, L = 3 £ , £ c = 5 £ in the absence of the field ( £ c = Hvf/2ttT, vf = HKp/M, £ = ^/h/mu>o ). 
This choice of W leads to N = 98, and \A(u m ,H)\ the pair amplitude in superconductors is 
assumed to be unchanged by a magnetic field. 

From this result, we see that, in a weak magnetic field oo c /ujq < 1, the Josephson current shows a 
dependence similar to a Fraunhofer pattern as Josin(7r$/<I > o)/(7i" < ^/ < ^o) with $ being the magnetic 
flux in the normal metal and a flux quantum <3?o = hc/2e, although the distance between successive 
maxima gradually increases as uj c /ujq exceeds about 0.5 . By fitting this numerical result with the 
above formula of a Fraunhofer pattern, we can estimate the effective width of the normal metal to 
be about 26 £, while the classically defined width W at H = is 28 £. 

These numerical results are understood as follows based on eq. (2.16). In such a weak field as 
LUc/^o < I, the integrand in eq. (2.16) has a strong peak around q = — (a\ + 02) with the width 
1 61 — 62 1- Other region of q does not contribute to the integrated value so much, because of the 
oscillation. In order to understand the dependence on parameters qualitatively, we approximately 
replace this integration by the value at q = — (a\ + 02) multiplied by | b\ — 62 1- 




(3.1) 
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(ai + a 2 )lo 
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(4.1) 



2biL - ft 2 | g 2 cos[(ai - a 2 )L] 

-F ni „ 2 [z] has the maximum around \z\ = l/2\/2iV • (M — Af + 1/2) • 7T7/2 where M = 
max{ni, 712}, -/V = min{ni,n2} and 7 = 0.6 (for details, see Appendix C). Therefore, in the sum- 
mation over ni and n 2 , contributions from terms satisfying the condition 



7T7 



27i (ai+a2) "27^r-r 2 - ni+ 2j 2 

1 n (a\- a 2 ,) 7T7 



2^/2^" ftw Kj, 2 ' 
dominate. For such pairs, the argument of the cosine in eq. (4.1) is given by 



(4.2) 



H-L-W cS ( ni ) 

(01 - a 2 )L = vr (4.3) 



where <l?o = /ic/2e and W e g(ni) is defined by 



WeflF(ni) = — .Wdl + 

7T7 V VcjoZ V A* 



W / cflf( n i) can be considered as an effective width of the channel n\. This effective width H / C flf( n i) i s 
an increasing function of ri\ and for the largest n\ ( huj{n\ + 1/2) ~ /j ), W e g(wi) approximately 
coincides with the classically estimated width W. From eq. (4.3), each maximum of the oscillation 
as a function of H corresponds to the change of the number of quantum flux <&o penetrating the 
region which electrons in the channel n\ effectively occupy. After taking the summation over m, 
we get a Fraunhofer pattern. 
In the absence of a magnetic field, because of the following relation 



lim F nin2 
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S ni n 2 , (4.5) 



the Cooper pairs consisting of electrons in the same channel as ( —kp(ni), fcp(rii) ) mainly con- 
tribute to Josephson current, and the Cooper pairs have a zero momentum as shown in Fig. 6(a). 
The above analysis shows that, as a magnetic field is applied in the region uj c /loq < 1, the pair 
states like (— £^(712), /cf(tji)) becomes effective where {ni,n 2 } meets the condition eq. (4.2), and 
the Cooper pair has a momentum, pF{n\) —pF{ n 2) (see Fig. 6(b)). Therefore, it can be concluded 
that as a magnetic field is applied the Cooper pairs consisting of electrons in different channels be- 
gin to contribute and the pairs begin to have a nonzero momentum. This change of the electronic 
states in the pairs contributing to the Josephson current result in a Fraunhofer pattern. 

The fact that the Cooper pairs contributing to the current begin to have a finite momentum leads 
to the spatial separation of Cooper pairs flowing in the opposite direction each other, although there 
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is a much spatial overlap under a weak magnetic field. This tendency can be easily understood 
by considering that the center of the single-particle wave function in the y-direction is given by 
yo(k) = Mq/Ic, proportional to the wave number k in the x-direction as discussed in §2. 

The increase of the difference in field between successive maxima at uj c /uo > 0.5 in Fig. 5 can be 
understood as follows. The difference between n\ and ri2 satisfying eq. (4.2) becomes larger as the 
field becomes stronger. However, the number of channels in the normal metal is restricted to about 
li/%uj. Therefore as the field is increased, the number of Cooper pairs contributing to Josephson 
current decreases and the contribution of the Cooper pairs with a small n remains. This leads to 
the increase of the distance between successive maxima as is seen by eq. (4.3). 

§5. Dependences in the Region of Strong Magnetic Field and Large N 

We now consider the effect of a strong magnetic field (uj c /ujq > 1). The number of channels are 
assumed to be large as in §4 (N = 98). 

The numerical result is shown in Fig. 7. All parameters are the same as those in §4 and the 
pair amplitude \A(u m , H)\ is assumed to be unchanged by the field. As the cyclotron frequency uj c 
increases beyond Uo, Josephson current starts to exhibit another type of oscillation which eventually 
dominates the Fraunhofer pattern in the strong field limit. The distance between successive maxima 
is not so different from that of a Fraunhofer pattern. 

Such behaviors can be understood as follows. In eq. (2.16) the peak at q = — (a\ + 02), which 
was important under the weak field, does not contribute in the present case of uj c /ujq >1, because 
F nin2 [qlQ / 21"^} at q = —(01 + 02) is so small. Instead, in the present case of co c /loo > 1, contributions 
from small z is important and we approximate F ni1l2 [z] as follows, 

{ Q[m] 2 / _ x 

[ |*| < QH ], (5.1) 
(n 1 ^n 2 ) 

where Q[n\\ = \/2yJ2n\ + 1 . Then we can perform the (/-integration in eq. (2.16) to obtain 



2 ai L - tan" 1 



e -(ci+26i)L 



T|A(^)p(f)S^^.£^, ,, 2 ) 



where c\ = 2/^//q • and we assumed the relation c\/2a\ <^ 1 and c\L < 1. 

This evaluation implies that Josephson current is mainly carried by the Cooper pairs consisting 
of electrons (kF(ni), (see Fig. 6(c)). This tendency can be understood by considering 

that under such a strong magnetic field the center of the single-particle wave function in the in- 
direction yo(kF(ni)) gets close to the classically estimated edge W/2, and the spatial extent in the 
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y-direction Iq becomes small. In this case, as seen from eq. (2.16), a main contribution comes from 
states satisfying n\ = n<i and yo(^i) — 2/0(^2)) namely k\ ~ £?2 since yo(k) oc k. Therefore Cooper 
pairs flow through the edge state in the normal metal as schematically shown in Fig. 8. 

Each term in eq. (5.2) involves an oscillatory function of H, cos(2aiL). This argument, 2a\L, 
can be expressed as follows, 

2 



o r H.L.2y Q ( ai )^ 

2 ai L = ir (5.3) 

<±>o 

where 2yo(ai) — 2yo(^F(^i)) = V§{kp(n\)) — yo(— fcp(ni)), i.e. the distance between two edge states 
around y = yo(kp(ni)) and y = yo(— fc_F(ni)), and the factor (uj/uj c ) 2 approaches to 1 in the large 
H limit. Therefore this argument corresponds to the number of the flux quantum $0 surrounded 
by the edge states as shown in Fig. 8. This oscillation results from the Aharonov-Bohm interference 
between the Cooper pairs, and this situation coincides with the result by Ma and Zyuzinljl under 
a strong magnetic field, the Josephson current is carried by Cooper pairs which consist of electrons 
having the same velocities and flow along the edge of the two-dimensional normal metal, which 
leads to the Aharonov-Bohm type oscillations as in quantum dots. 

In this system with many channels, as a magnetic field is increased, Cooper pairs rearranges 
themselves from (— kp(n2), A:_p(ni)) to (fcp(ni), fei?(ni)) in due order from Cooper pairs with a larger 
n. This is why Josephson current starts to show another type of oscillation around uj c /ujq ~ 1. 
In the strong magnetic field limit, almost all the Cooper pairs finally flow through the edge states. 
This leads to a complicated oscillation because the amplitude and the periodicity of the oscillatory 
behavior in each channel are not correlated as is seen in eq. (5.2). However the distance between 
successive maxima of Josephson current are not so different between weak and strong field limits, 
because the argument in eq. (5.3) is approximately the same form as one in eq. (4.3), although 
under a high magnetic field the chemical potential /x decreases (i.e. the width W decreases). 

§6. Quantum Limit ( Small N ) 

In §4 and §5, we assumed that there are many channels (N ~ 100) in the normal metal, i.e. 
the width W is large enough. In this section, we consider the case where there are not so many 
channels in the normal metal. 

The numerical results for several choices of W are shown in Fig. 9. In order to see the dependence 
on N clearly, the width W in Fig. 9(a) is chosen as 20. 1£ corresponding to N = 51 which can be 
considered in the middle of a large N and a small N. Here, \A(uj m , H)\ is assumed to be unchanged 
by the field as before. As N decreases, we see an appreciable deviation from a Fraunhofer pattern 
in the weak field region (uj c < ujq) and a less oscillatory behavior due to the edge states in the 
strong field region (u c > ujq). 

In order to understand these properties, we consider a one-channel system (N = 1) as a limiting 
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case, where the Josephson current is given 



e-* L (ft)' + 



COS 


2a,L - tan- 1 (ft) 
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'(a? + b\) ■ 


f i+ te) 


2 ) 





3 -(ci+26i)L 



Cl 



\2ai y 

(6.1) 

where ci = 2l 2 c jl\ ■ Q[ni]. Here we approximated F„ ini [z] in eq. (2.16) by the Lorentzian as in eq. 
(5.1). 

Since the factor c\/2a\ is larger than 1 in the region uj c /ujq < 1 (ci/2ai = oo when oj c — ► 0), the 
first term dominates, which means that Josephson current does not show any oscillatory dependence 
in the region of the weak field. As shown in Fig. 6(b), the Fraunhofer pattern is caused by Cooper 
pairs consisting of electrons belonging to the different channels. Therefore as TV decreases, the 
interference between Cooper pairs is weakened. 

Under the strong field {uj c /lvq > 1), the factor c\/2a\ becomes smaller than 1, and the 1st term 
becomes small. However the factor c\L in the damping factor exp(— (c\ + 2b\)L) in the second term 
is large for a small n\, so the second term is much smaller than the first term and the oscillatory 
behavior does not appear. 

Especially in Fig. 9(b) and (c), the decrease of the number of channels as the field increases 
affects Josephson current. This is because the contribution of each channel to Josephson current 
becomes large in such a narrow system. 

It can be concluded that as the number of channels decreases the field dependence in the weak 
field region (lo c /ujq < 1) deviates from a Fraunhofer pattern and the amplitude of the oscillatory 
behavior due to the edge states under the strong field (uj c /ojq > 1) is reduced. 

§7. Summary and Discussion 

We have studied the effect of an external magnetic field on Josephson current in a SNS system 
where two bulk superconductors are attached to a confined two-dimensional normal metal and the 
temperature is assumed to satisfy £ c < L (£ c is the coherence length of the pair in the normal 
metal). We first studied the case that the normal metal has many channels (N ~ 100). In the 
absence of field, Cooper pairs which consist of two electrons in the same channel and have an 
essentially zero momentum contribute to the current. For the weak field (oj c /ujo < 1), however, 
the contribution of pairs which consist of electrons in the different channels and have the finite 
momentum pF{n\) — pf^) becomes larger. This leads to a Fraunhofer pattern. As the field is 
increased further {uj c /ujq > 1), Cooper pairs with a larger n rearrange themselves into the states 
that two electrons are in the same channel and have the same momentum. This means that 
Josephson current flows through the edge states, and shows an oscillatory dependence due to the 
Aharonov-Bohm interference between Cooper pairs. This is physically the same as the result by 
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Ma and ZyuzinE|)0 who studied the case where two superconducting leads are attached via point 
contacts to a confined two-dimensional normal metal and the edge states of the lowest Landau level 
are filled. 

We next examined the dependence on N the number of channels in the normal metal at H = 0. 
The decrease of N weakens an interference between Cooper pairs even in the region of the weak 
field and leads to a deviation from a Fraunhofer pattern. In a limiting case, Josephson current 
in an one-channel system (JV = 1) does not show a Fraunhofer pattern dependence at all. Under 
the strong field, the amplitude of the oscillatory behavior due to the edge states is reduced as N 
decreases. 

In this paper, we neglected Zeeman energy which splits a Fermi wave number of the each channels 
as kp + Ak, kp — Ak. Under the strong field, the single-particle wave function ipnk{ r ) is localized 
around the center yo(k) in the y-direction with the spatial extent l c y/n + 1/2. In the absence of 
the Zeeman energy, a much contribution comes from the pair states where two electrons belong 
to the same channel and the centers of two single particle wave functions in the y-direction are in 
an almost same position. The Zeeman energy separates this overlapping two single-particle wave 
functions in the y-direction as yo(kp + A) and yo(kp — A). Therefore Zeeman energy will reduce 
Josephson current. 

In this paper, we assumed that the confining potential caused by the gate voltage is parabolic. 
The self-consistent calculation for a GaAs/AlGaAs heterostructure by Laux et aB sug gests that 
the parabolic potential is appropriate for a system with one or two channels (N = 1 or N = 2), but 
as charge accumulates in the well (i.e. N is increased in the region ), the potential gets close to 
a square well type. If this is the case, the behavior of Josephson current in the system with a few 
channels as discussed in §6, can be considered rather realistic. In the system with a large iV(~ 100) 
as discussed in §4 and §5, however, the edge states may be located nearer to the classically estimated 
edge ±W/2 regardless of n. Therefore the periodicities of the oscillations due to edge states in eq. 
(5.3) can be almost the same for all channels, and the oscillatory dependence due to edge states 
is expected to become simpler than one shown in Fig. 7. On the other hand, Fraunhofer pattern 
under the weak field is considered not to be affected so much by this alternation of a confining 
potential. From these discussions, it can be expected that the field dependence in a strong field 
region will reflect the actual shape of the confining potential. 
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Appendix A: Symmetries of the Green's Function 

In the derivation eq. (2.15), we have used the following symmetries of the Green's functions; 

(A-l) 



GNa(ri,r r ; -u m ) = GNa{r r , ri;u m )* 

dyi J dy r G N ^{r u r r ;LOm)Q Ni (r u r r ]-u m ) 

dyi / dy r Q N ^(r r ,rr,uj m )Q N i(r r ,ri;-uj m ) 



(A-2) 



Appendix B: Calculation of the Green's Function 

By use of the formula,© 

/oo „ , 

dx e~ x H m [x + y]H n [x + z] = 2 n ^ m \z n ~ m L^- m) [-2yz] 
-oo 

we can transform J c as the following, 

Jc = j^-T |A(w m ,#)| 2 Re J dyi J dy r G N ^(r u r r ;u m )g Nl (r u r r 



(n>m), (B-l) 



-0J r . 



u; m >0 

dK 



211 



AKL 



(B-2) 



2vr (K + q - 2n l ){K + q + 2k x ){K -q- 2k 2 ){K -q + 2k 2 ) ' 

where K = k\ + k 2 (the total momentum) , q = k\ — k 2 (the relative momentum) , and k\ and k 2 are 
treated as continuous here by the following reason. At the interfaces between the normal metal and 
the superconductor, there will be a repulsive potential, i.e. potential barrier, in the x-direction, 
through which electrons in the normal metal region are coupled to the propagating states in the 
bulk superconductors. After performing the K-integration, we obtain eq. (2.16). 

Appendix C: Estimation of Maximum Point of F niri2 [z] 

The asymptotic behavior of the Laguerre polynomial© is given by 

1\ vr 



1 i£ _ £L 1. £i 1. 

e 2 x 2 4 n2 4 cos 

7T 



2Jnx 



a + 



2 2 



+ 0(n 7 < 



(C-1) 



which is adequate for a large n and roughly 

l/2y/n • (a + 1/2) • vr/2 < ^ < 1/2^/n • (a + 1/2) • vr/2 + ^tt/2. 
By use of this, F nin2 [z] is approximated by 



: COS 



2VN ■ 2z 2 -IM-N + - 



(C-2) 



ttVN ■ 2z 2 

where M = maxjni, n 2 }, N = min{ni,7i2}. 

From this asymptotic form, F niH2 [z] can be expected to have a maximum around \z\ = l/2y/2N ■ 
(M — N + 1/2) • 7r/2. However, this asymptotic form is valid in the larger region than this maximum 
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point. Therefore the real maximum point deviates from the expected one and can be written as 



where 7 is weakly dependent on n l5 7J2 and ranges roughly from 0.5 to 0.7. For simplicity, we have 
used the average value of 0.6 for 7. 
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Fig. 1. Schematic illustration of SNS system considered in this paper. The origin of the coordinate system (x,y) is 
chosen at the center of the left SN interface. 



Fig. 2. Dispersion relation, E n k vs k for electronic subbands arising from electronic confinement in zero magnetic 



Fig. 3. The lowest order process of Cooper pair propagation considered in this paper. Electronic state in the normal 
metal is specified by {n, k}. 




(C-3) 



field. 
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Fig. 4. Schematic picture showing the dependence of characteristic magnetic field effects on N the number of 
channels in the absence of the field. 



Fig. 5. The dependence of Josephson current on a weak magnetic field (u> c /^o < 1) in a system with large N 
(N=98). The parameters are chosen as W = 28 £, L = 3 £, £ c = 5 £ (£ = yJU/muio). J c is normalized by 8emT/h 3 n- 
\A{ttT,H)\ 2 = Jo. 



Fig. 6. Cooper pairs giving a dominant contribution to Josephson current, (a)zero magnetic field, (b)weak magnetic 
field, (c)strong magnetic field. 



Fig. 7. Josephson current as a function of magnetic field in the case of strong magnetic field (w c /wo > 1). All 
parameters are the same as those in Fig. 5, i.e. W = 28 £, L — 3 £, £ c = 5 £ (£ = y/h/mujo) 



Fig. 8. Schematic illustration of the motion of Cooper pairs flowing through the edge state under a strong magnetic 
field. 



Fig. 9. Josephson current as a function of magnetic field for various choice of system width, W. (a) W = 20.1 £ 
(N = 51), (b) W = 9.2 £ (N = 11), (c) W = 5.3 £ (AT = 4), (d) W = 2.2 £ (AT = 1), where L = 3 £, £ c = 5 £, £ = 
\J%JmuJo~. 
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